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Abstract: By taking the high-energy limit of the two-loop and the three- loop four-point 
amplitudes in the maximally supersymmetric N = 4 Yang- Mill theory (MSYM), we test 
the validity of the loop expansion of the high-energy amplitude, beyond the next-to-leading- 
logarithmic (NLL) accuracy. We compute the three-loop Regge trajectory, and the two-loop 
and three-loop coefficient functions. These quantities are relevant for the BFKL evolution 
beyond NLL, as well as building MSYM two-loop and three-loop amplitudes with many legs 
in the high-energy limit, which in turn may be used as a powerful check of the evaluation 
of the corresponding exact amplitudes. 
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1. Introduction 



In the high-energy limit (HEL) s 3> \t\, any scattering process is dominated by the ex- 
change of the highest-spin particle in the crossed channel. Thus, in perturbative QCD the 
leading contribution in powers of s/t to any scattering process comes from gluon exchange 
in the t channel. Building upon this fact, the BFKL theory models strong- interaction 
processes with two large and disparate scales, by resumming the radiative corrections 
to parton-parton scattering. This is achieved to leading logarithmic (LL) accuracy, in 
ln(s/|i|), through the BFKL equation [1, 2, 3], i.e. a two-dimensional integral equation 
which describes the evolution of the t-channel gluon propagator in transverse momentum 
space and moment space. The integral equation is obtained by computing the one-loop LL 
corrections to the gluon exchange in the t channel. They are formed by a real correction, 
the emission of a gluon along the ladder [4, 5, 6], and the leading virtual contribution of 
the gluon loop, the one-loop Regge trajectory [1]. The next-to-leading-logarithmic (NLL) 
corrections to the BFKL equation have been computed as well [7, 8]. The virtual part of 
the NLL kernel is provided by the two-loop trajectory [9, 10, 11, 12, 13], whose evaluation 
is based upon the assumption of Regge factorization of the scattering amplitudes beyond 
leading logarithmic accuracy. 

Recently, Bern, Dixon and Smirnov (BDS) have proposed an ansatz [14] for the £-loop 
n-gluon scattering amplitude in the maximally super symmetric N = 4 Yang-Mill theory 
(MSYM), with the maximally- helicity violating (MHV) configuration and for arbitrary I 
and n. In Ref. [14], an analytic expression for the exact three-loop four-gluon amplitude, 
as well as iterative relations for the n-point MHV MSYM amplitudes, have also been 
provided. Using the results of Ref. [14], we test the high-energy factorization of the gluon- 
gluon scattering process. The paper is organised as follows: in Section 2, we analyse 
high-energy factorization at three-loop accuracy in perturbative QCD; in Section 3, we 
take the exact two-loop and three-loop four-point MSYM amplitudes in the HEL, and 
derive the three-loop Regge trajectory, as well as the two-loop and three-loop coefficient 
functions; in Section 4, we make contact between these quantities, the BDS ansatz and the 
iterative relations for the MHV MSYM amplitudes; in Section 5, we sketch an outlook for 
the BFKL evolution in MSYM beyond NLL accuracy. 



2. High-energy factorization 



In perturbative QCD, the simplest process is parton-parton scattering, which in the HEL 
occurs through gluon exchange in the t channel. It is convenient to focus on gluon-gluon 
scattering. In the HEL, the tree-level amplitude for g a gb — ► 9 a' 9v ma Y De written as [1], 



igsf aca 'c^(p a ,p a ,)] \ ig s f bcb 'c^( Pb ,p b/ ) 



(2.1) 



where a, a', b, b' represent the colours of the scattering gluons. The gluon coefficient func- 
tions C^°\ which yield the LO gluon impact factors, are given in Ref. [1] in terms of their 
spin structure and in Ref. [15, 16] at fixed helicities of the external gluons. 
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The colour decomposition of the tree-level n-gluon amplitude in a helicity basis is [17] 



M (0) =2 n/2 g n-2 £ t r(I<M • ■ • 2<(»> ) (p ff(1) , Ml) ; P<T(n) , v a{n) ) , (2.2) 

where di,...,d n , and i/i,...,i/ n are respectively the colours and the polarizations of the 
gluons, the T's are the colour matrices* in the fundamental representation of SXJ(N) and 
the sum is over the noncyclic permutations S n /Z n of the set [l,...,n]. We take all the 
momenta as outgoing, and consider the MHV configurations ( — ,—,+,...,+) for which the 
tree-level gauge-invariant colour-stripped sub-amplitudes, rrin\pi, v\\ ...;p n , v n ), assume 
the form 

m£H- ~ +■ ■■■> +) = 7 v w r ■ (2-3) 

(P1P2) ■ ■ ■ (Pn-lPn)(PnPl) 

where i and j are the gluons of negative helicity. In gluon-gluon scattering, only the helicity 
configurations (—,—,+,+) occur at tree level, and out of the six possible colour configu- 
rations only four are leading in the HEL [18]. They are those corresponding to s-channel 
helicity conservation. For instance, let us label the gluons clockwise and consider the he- 
licity configuration (b—,a—, a'+, b'+), with a and b incoming and a' and b' outgoing, then 
the four leading colour configurations are (b,a,a',b'), (b,b',a',a), (b,b',a,a'), (b,a',a,b'). 
The latter two, corresponding to the helicity ordering (—,+,—,+), can be obtained from 
the former two, corresponding to the helicity ordering (— , — , +, +), by s «-> u channel ex- 
change. In addition, in the HEL rn^\— , +, — , +) = — rn^\— , — , +, +), thus the different 
colour configurations contribute to Eq. (2.2) with alternating signs, in such a way that the 
traces of T matrices combine to form the structure constants of Eq. (2.1) [15]. 

The virtual radiative corrections to eq. (2.1) in LL approximation are obtained, to all 
orders in a s , by replacing [1] 



(2.4) 



in eq. (2.1), where a(t) can be written in dimensional regularization in d = 4— 2e dimensions 
as 

a(t)=g 2 s cr (V) N ~> ( 2 - 5 ) 



-t 



with N colours, and 



i r(i + e )r 2 (i- £ ) 
cr = Wf^ r(i - 26) • (2 ' 6) 



The fact that higher order corrections to gluon exchange in the t channel can be accounted 
for by dressing the gluon propagator with the exponential of Eq. (2.4) is called the gluon 
reggeization, with a(t) the Regge trajectory. In order to go beyond the LL approximation, 
we need a prescription that disentangles the virtual corrections to the coefficient functions in 
Eq. (2.1) from those that reggeize the gluon (2.4). The prescription for doing so is supplied 
by the general form of the high-energy amplitude for gluon-gluon scattering which arises 



*We use the normalization tr(T c T d ) = S cd /2, although it is immaterial in what follows. 
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from a single reggeized gluon exchanged in the crossed channel [19] 



M 4 



ig s f aca C(p a ,p a , 



—s 



<*{t) / s \ a(t) 
+ i ~t 



ig s f bcb 'C( Pb , Pb , 



(2.7) 



The first Regge trajectory in Eq. (2.7) corresponds to the s-channel physical region; the 
second to the u channel. The s and u channels have different analytic properties, how- 
ever Eq. (2.7) presumes that the corresponding sub-amplitudes still differ only by a sign 
m 4 (— , +, — , +) = — to 4 (— , — , +, +). It is clear that the equals sign in Eq. (2.7) cannot be 
expected to hold strictly, but, in fact, holds only up to NLL accuracy [20]. However, we 
show that, at the colour-stripped amplitude level, the following high-energy prescription is 
valid beyond NLL accuracy, 



m 4 (-, -,+,+) =m% = s[g s C(p a ,p a >)] - — [g s C(p b ,p b >)] , 



a(t) 



in the s-channel physical region, and similarly 



m 4 (-, +,-,+) =m\ = s[g s C(p a ,p a ,)] - ( — 



a(t) 



[gsC{p b ,p b < 



(2.8) 



(2.9) 



in the u-channel physical region. In Eqs. (2.7)-(2.9), the gluon Regge trajectory has the 
perturbative expansion, 



fj_ 
-t 



2 \ e 



(2.10) 
(2.11) 



a(t) = ~g 2 s (t)aW + ^ (t)a (2) + ^6 (t)a (3) + 0( ~8 } 
with the rescaled coupling 

g 2 s (t) = gW 

and with = 2N/e given in Eq. (2.5). The coefficient functions C can be written as 

C = C(°>(1 + ~g 2 s {t)C^ + t s {t)C {2) + f s {t)C {Z) ) + 0(f s ) . (2.12) 

Because the coefficient function C is real (up to overall complex phases in induced 
by the complex-valued helicity bases), Eqs. (2.7)-(2.9) imply that the imaginary part of 
the amplitude comes entirely from the Regge trajectory in the s channel, according to the 
usual prescription ln( — s) = ln(s) — m, for s > 0. 

The expansion of Eqs. (2.8) and (2.9) can be written as, 



r6 m i(3) 



m 4 = ml I 1 + g s + g s ml + g s ml + 
with i = s,u. The one-loop coefficients of Eq. (2.13) are, 

m ; (1) = a«L+ 2CW, 



(2.13) 



m 



s(l) tt(l) • (1) 



(2.14) 



with L = ln^^^. The one- loop trajectory, qW, is universal, i.e. it is independent of 
the type of parton undergoing the high-energy scattering process. It is also independent 
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of the infrared (IR) regularisation scheme. Conversely, the one-loop coefficient function, 
C(i) , is process and IR-scheme dependent. was computed in conventional dimensional 
regularization (CDR)/'t-Hooft-Veltman (HV) schemes in Ref. [19, 20, 21, 22, 23], and in 
the dimensional reduction scheme (DRED) in Ref. [20, 23]. 
The two-loop coefficients of Eq. (2.13) are, 



m 



«< 2 > = I ( a m 



L 



+ 



(«( 2 >+2C( 1 )«( 1 )) L 



(2.15) 



+ 2C^ + (C^ 



m 



s(2) 



u{2) K ( (l)V 



— ITT 



(««) 2 L + «( 2 )+2CW a W 



(2.16) 



The two-loop trajectory, oft\ was computed in the CDR scheme in Ref. [9, 10, 11, 12, 13]. 
We note that the coefficients of the double and the real parts of the single logarithms are 
the same in the s and the u channel. Therefore, it is correct to use Eq. (2.7) at NLL 
accuracy. However, the real part of the constant term is channel dependent, and thus care 
must be used in extracting the two-loop coefficient function from it^. 
The three-loop coefficients of Eq. (2.13) are, 



m 



u(3) 



+ 



P (aW+cWaW)L 2 

a (3) + 2a (2) c (i ) + a (i) f 2 c<n + (cw) 



(2.17) 



+ 2C (3) + 2C (2) C (1) , 



m 



s(3) 



m 



u(3) 



( (l)^ 3 



— ITT 



(2.18) 



We note that the coefficients of the triple and the real part of the double logarithms are 
the same in both s and u channels, thus maintaining the correctness of Eq. (2.7) at NLL 
accuracy. However, although the three-loop trajectory is universal, the coefficient of the 
single logarithm, and thus the extraction of the three-loop trajectory, depends on the 
channel under consideration. 



tin Ref. [13], the constant terms for the gluon-gluon, quark-quark and gluon-quark scattering processes 
were evaluated by projecting the two-loop amplitude on the tree amplitude, which entails averaging over 
the s and u channels. However, an apparent discrepancy in the factorization between the quark and the 
gluon amplitudes prevented any conclusions from being drawn about the two-loop coefficient functions. 
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3. The MSYM amplitudes 



The QCD four-point parton-parton scattering amplitudes are known at two- loop accu- 
racy [24, 25, 26, 27, 28, 29, 30, 31, 32]. In QCD, the gluon loop can be decomposed 
into a MSYM multiplet, an N = 1 chiral multiplet, and a complex scalar. The MSYM 
contribution is the simplest to evaluate, and captures most of the IR behaviour of the full 
gluon loop. In fact within MSYM, the gluon-gluon scattering amplitude is known at three 
loops analytically [14] and at four loops through 0(1/ e) numerically [33]. Bern, Dixon and 
Smirnov have proposed an ansatz [14] for the MHV /-loop n-gluon scattering amplitude for 
arbitrary I and n. The ansatz agrees with the direct evaluation of the three-loop four-point 
amplitude* In addition, in Ref. [40, 41] it has been shown that in the HEL the BDS ansatz 
for the four-point amplitude exhibits the Regge behaviour of Eqs. (2.8) and (2.9). Thus, 
we shall use the direct evaluation of the three-loop amplitude and the BDS ansatz to derive 
the relevant quantities for the three-loop high-energy factorization in MSYM. 

We start with the HEL of the QCD one-loop colour-stripped amplitude for gluon-gluon 
scattering, which is known to all orders in e [23], and use the maximal trascendentality 
principle [42] to select the MSYM contribution. Using the conventions of Eq. (2.13) and 
rescaling the coupling (2.11) as 

f s {t) = ~g 2 s (t)N , (3.1) 



we obtain 



4 2L 2 n , T 4 2 o, 3 2tt 6 

■? + T + 7r +2C3£ + T5 e +2C5£ + 315 ( 



yYlf*^^ == TYl^^~^ (3 3^ 

4msym ^msym e V • / 

From Eq. (3.2), we see that the MSYM one-loop trajectory is precisely the same as that in 
QCD, Eq. (2.5), while the MSYM one-loop coefficient function is, to all orders in e, given 
by 

^(i) _ ^(l + e)-2V>(-6)+^(l) 
msym — ■ y^-V 

At two loops, we take the direct calculation of the exact MSYM amplitude, given in 
Ref. [14] to 0(e 2 ), and evaluate the HEL in the u and s physical channels, 

u(2) 8 8L / 2 Hvr 2 \ 1 (5tt 2 \1 



'"The BDS ansatz was first predicted to fail by Alday and Maldacena [34], for amplitudes with a large 
number of gluons in the strong-coupling limit. They claimed that the finite pieces of the two-loop six gluon 
amplitude would be incorrectly determined. This prediction was backed up by Drummond et al. [35], who 
considered the finite contribution in the dual theory [36] by computing the hexagonal light-like Wilson 
loop at two loops. The conclusion was that either the BDS ansatz is wrong, or the equivalence between 
Wilson loops and scattering amplitudes does not work at two loops. Recent numerical results for the finite 
part of the MHV six-gluon amplitude in MSYM [37] have confirmed the equivalence with the finite part 
of the light-like hexagon Wilson loop [38] thereby disproving the BDS ansatz. Furthermore, the analytic 
structure of the two-loop six-gluon amplitude in the multi-Regge kinematics is expected to be at odds with 
the structure of the BDS ansatz [39]. 
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1 ^ L -Vo" 



/2ir 4 



4ir 2 



;, + ll^ L + - C3 - 86C5)e 



(3.5) 



+ 



(6vr 2 C 3 + 86C5) L - 94(1 - 



777r b 



m 



8(2) 

4msi 



u(2) 8i7T 
m 4:MSYM + e 3 



180 

2ir 2 + 4mL 5in 3 



e 2 + 0(e 3 ) , 
2C 3 ivr 



3c 



2ivr 



e - (6C 3 vr 2 + 86C5) «vre 2 + C(e 3 ) . 



u(2) 
^MSYM 



and m 



«(2) 

&MSYM 



(3.6) 



differ by 



As expected from Eqs. (2.15) and (2.16), the real parts of m 

the constant term 2ir 2 /e 2 . Note that the average over the s and u channels is in agreement, 
to 0(e°), with the terms of highest transcendentality of the HEL of the projection of the 
two- loop amplitude on the tree amplitude [13]. 

Equating the coefficients of the single logarithm of Eqs. (2.15) and (3.5), we obtain the 
MSYM two-loop trajectory [42] in the DRED scheme, 



= - 2C 3 - + (6A3 + 82C 5 )6 2 + 0(e 3 



(3.7) 



which is in agreement, to O(e ), with the terms of highest transcendentality of the QCD 
two-loop trajectory [13]. 

Equating the coefficients of the constant term of Eqs. (2.15) and (3.5), we obtain the 
MSYM two-loop coefficient function, 



C 



(2) 

MSYM 



57T 2 1 C3 

e 



- - -K 1 



72 

95 , 2 113vr 6 \ 
2 



e 2 + C(e 3 ) . 



(3.8) 



The evaluation of Eq. (3.8) has been performed also in the s channel through Eqs. (2.16) 
and (3.6), thus checking the consistency of Eqs. (2.8) and (2.9) at two loops. 

At three loops, we take the direct calculation of the exact MSYM amplitude, provided 
in Ref. [14] to 0(e°), and evaluate the HEL in the u and s physical channels. We find 



m 



u(3) 

^MSYM 



32 16L 
+ — =- 



- 8L 2 



207T 2 



+ 



3e 6 



181 



+ -L 3 - 6vr 2 L + 8C3 



26 , 112 



1 

952 



^ - — ttU + — C 3 vr z - — Cs - 



- ^rL 2 -i 1.8J V) + 



45 



196tT 2 

9 



135 
3284 



c! + 



27 
88747tt 



6 



51030 



+ 000, 



(3.9) 



m 



*(3) 



'MSTM 



m 



u(3) 



*MSYM 



Wiir 8vr 2 + 16Livr 
5 I ^ 



4vr 2 L + 4L 2 - 



22vr z 



ITT 



4ir 4 



+ 



2tt 6 

+ -T + 

45 



87T 2 

3 

4tt 4 



L - 4C 3 ivr 



1 26 iir 5 

7 + 135 ~r 



)5 L+^C 3 7r 2 + 484C 5 ) tt/ f £){<). 



(3.10) 
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The coefficients of the double logarithm of Eqs. (2.17) and (3.9) provide a cross check of 
the two-loop trajectory a^ 2 \ while comparing the coefficients of the single logarithm we 
obtain the MSYM three-loop trajectory, 

a (3) = §^ + ^C3 + 16C 5 + 0( e ), (3.11) 

in agreement with Ref. [39, 40]. The trajectory was also implicitly evaluated in Ref. [41]. 
Equating the coefficients of the constant term of Eqs. (2.17) and (3.9) we find the MSYM 
three-loop coefficient function, 



( 3 ) _ 4 2vr 2 1 217vr 4 1 / llvr 2 224 \ 1 
msym ~ ~3^6 + ~3~^4 + ~sw"^ + \ 27~ + ~3~ ) 1 



+ 



f 796 211861^ 
I 9 S3 408240 J 



The evaluation of Eq. (3.12) has been performed also in the s channel through Eqs. (2.18) 
and (3.10), thus checking the consistency of Eqs. (2.8) and (2.9) at three loops. 

4. The Bern-Dixon-Smirnov ansatz 

The BDS ansatz prescribes that the n-point MHV amplitude be written as, 



mffl exp 



L=l 

oo 

a 1 (/«(e)M«(Ze) + Const® + E®(e)) 



i=i 



(4.1) 



(4.2) 



where 

a = ^Le~^ 
(4vr)2- e 

is the 't-Hooft gauge coupling, and with 

f {l) (e) = f® + ef® + e>f®, (4-3) 

where f^\e) = 1, and f® is proportional to the Hoop cusp anomalous dimension, 7^ = 
and f® is related to another quantity, Q® = 2f® /I, which enters the IR Sudakov 
form factor and accounts for virtual divergences which are not simultaneously soft and 
collinear [43, 44]. In Eq. (4.1), Const® are constants, and E®(e) are 0(e) contributions, 
with Const^ = and En\e) = 0, and M^\e) is the L-loop colour-stripped amplitude 
rescaled by the tree amplitude. In the convention and notation of Eq. (2.13), the four-point 
amplitude is given by, 
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with 



G(e) = 



e -ie _ 2e ) 



l + 0(e 2 



(4.5) 



r(l + e)r 2 (l-e) 

Thus, using the rescaled coupling (3.1), the BDS ansatz (4.1) for the four-point amplitude 
becomes 



ni4 



m 



(0) 



i + ^f(i)4 L) ( e ) 



(0) 

= m\ exp 



L=l 

oo 



(2G(e))' + Const® + 0(< 



(4.6) 



Eq. (4.6) applies equally well in either the s or the u channels. Substituting the HEL 
one-loop amplitude (2.14) in Eq. (4.6) and comparing with the expansion (2.13) of the 
high-energy factorization (2.8) or (2.9), we see that the coefficient of the single logarithm 
allows us to read off the value of the gluon trajectory, 



a 



(i) 



- e G(e) 



a 



(2) 



(4.7) 



a 



(3) = - 



3e 



and in general 



aU = Z -G\e)fU(e) 



(4- 



From Eq. (4.8), we see that only the first two terms of the /®(e) function (4.3) enter the 
evaluation of the Regge trajectory. Dropping G(e), which does not contribute in Eq. (4.7) 
to O(e ), and using the and /® functions [14], 



/ (2) (e) = -C2 " Cse - Ue 2 , 

/ (3) (e) = y U + (6Cs + 5C 2 C 3 )e + (ciCe + c 2 d)e 2 , 



(4.9) 



we see that Eq. (4.7) agrees with Eqs. (3.7) and (3.11) to 0(e°). The coefficients ci , C2 are 
unknown, but they do not enter the evaluation of the Regge trajectory. 

Through the iterative structure of the MSYM amplitudes, it is possible to express 
the coefficient functions at a given loop in terms of the coefficient functions at a lower 
number of loops. Using Eq. (4.6), the iterative structure of the two-loop four-point MSYM 
amplitude [45, 46] is given by 



1 

2 L 



2 2G 2 (e) 
+ ~G{2e) 



( e ) m W (2e) + 4 Const( 2 ) + O(e) , (4.10) 



with Const^ = — Cf/2, and where the one-loop amplitude must be known to 0(e 2 ). Using 
the two-loop factorization, in either the u (2.15) or s channels (2.16), we find 



C 



(2) 



l r 



C { m ) sym^)\ + ^^/ (2) ( e ) G^sYM^e) +2ConstW + {e) , (4.11) 
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where the one-loop coefficient function is needed to 0(e 2 ). Eq. (4.11) agrees with Eq. (3.8) 
to O(e ). 

Similarly, the iterative structure of the three-loop four-point MSYM amplitude is [14], 
mf(e) = -\ \mf\e)] 3 + mf (e) (e) + -^^/( 3 )( e ) m<£\3e) +8 Const® + 0(e) , 



G(3e) 

(4.12) 



where m^\e) and m^\e) must be known to 0(e 4 ) and 0{e 2 ), respectively, and with 

^ <3, =(i + H c « + (-^ + H c| - (413) 

Using the three-loop factorization, in either the u (2.17) or s channels (2.18), we obtain 
the three-loop coefficient function 



clisYhti 6 ) — 3 C ( m ' S ym( € ) + G\J SYM {e)C y M > SYM {e) 

+ ^^/ (3) ( £ ) CmsymW + 4<W 3 > + O(e) . (4.14) 

The coefficients c±,C2 cancel when Eqs. (4.9) and (4.13) are used in Eqs. (4.12) and (4.14). 
Using the two-loop coefficient function to 0(e 2 ) (3.8), and the one-loop coefficient function 
to £>(e 4 ) (3.4), we see that Eq. (4.14) is in agreement with Eq. (3.12) to C(e°). 



1 



5. Conclusions 

The iterative structure and the exponentiated form of the MSYM amplitudes in the MHV 
configuration offer a useful computational lab to test high-energy factorization in the 
MSYM, and allow the derivation of some of the relevant quantities in the high energy 
limit. Using the MSYM two- and three-loop four-point amplitudes, we have tested the 
high-energy factorization of the colour-stripped amplitude. In particular, we have shown 
that it is valid beyond NLL accuracy, and we have verified that the factorization formulae 
(2.8) and (2.9) hold at three-loop accuracy. Accordingly, we have derived the three-loop 
Regge trajectory (3.11), as well as the two-loop (3.8) and three-loop coefficient functions 
(3.12). 

The three- loop Regge trajectory is one of the building blocks for deriving a BFKL 
evolution equation at next-to-next-to-leading logarithmic (NNLL) accuracy. The others 
are the tree vertex for the emission of three gluons along the ladder [47, 48], the two-loop 
vertex for the emission of a gluon along the ladder, and the one-loop vertex for the emission 
of two gluons along the ladder. In addition, if one wants to compute jet cross sections at 
a matching accuracy, the NNLO impact factors must be determined. The ingredients for 
that are the two-loop coefficient function evaluated here, as well as the tree coefficient 
function for the emission of three gluons [47, 48], and the one-loop coefficient function for 
the emission of two gluons. 

Furthermore the two-loop and three-loop Regge trajectories and coefficient functions, 
together with the vertices for the emission of one or more gluons along the ladder which 
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were not analysed in this work, can be used to build MSYM two-loop and three-loop 
amplitudes with a larger number of legs in the high-energy limit. This may serve as a 
powerful check on the structure of high multiplicity MSYM amplitudes. 
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